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Abstract 

The relativistic field theory model of the deuteron (RFMD) is re- 
formulated from the first principles of QCD. The deuteron appears as 
a neutron-proton collective excitation, i.e. a Cooper np-pair, induced 
by a phenomenological local four-nucleon interaction in the nuclear 
phase of QCD. The RFMD describes the deuteron coupled to hadrons 
through one-nucleon loop exchanges providing a minimal transfer of 
nucleon flavours from initial to final nuclear states and accounting for 
contributions of nucleon-loop anomalies which are completely deter- 
mined by one-nucleon loop diagrams. The dominance of contributions 
of nucleon-loop anomalies to effective Lagrangians of low-energy nu- 
clear interactions is justified in the large Nc expansion, where Nq is 
the number of quark colours. 
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1 Introduction 



In our recent publications [1-5] we have formulated the relativistic field 
theory model of the deuteron (RFMD) [1,2]. The first version of the RFMD 
[1,2] formulated in analogy with a phenomenological a-model had no obvious 
connection with QCD [2] . Unlike the potential model approach (PMA) [6-8] 
and the Effective Field Theory (EFT) approach [9-14] the RFMD takes into 
account non-perturbative contributions of high-energy (or short-distance) 
fiuctuations of virtual nucleon [N) and anti-nucleon {N) fields, i.e. NN 
fiuctuations, in the form of anomalies of one-nucleon loop diagrams^. The 
description of one-nucleon loop anomalies caused by contributions of high- 
energy (or short-distance) fiuctuations of virtual nucleon fields goes beyond 
the scope of both the PMA and the EFT approach. This is due to the 
absence in these approaches anti-nucleon degrees of freedom related to the 
non-perturbative quantum vacuum - the nucleon Dirac sea [15]. 

It is known that the nucleon Dirac sea cannot be ignored fully in the 
low-energy nuclear physics. For example, high-energy NN fiuctuations of 
the nucleon Dirac sea polarized by the nuclear medium decrease the scalar 
nuclear density in the nuclear interior of finite nuclei by 15% [16]. This 
effect has been obtained within quantum field theoretic approaches in terms 
of one-nucleon loop exchanges. Unfortunately, contributions of nucleon-loop 
anomalies have not been taken into account in these approaches. The RFMD 
allows to fill this blank. In fact, in accord the analysis carried out in Refs. [17] 
nucleon-loop anomalies can be interpreted as non-trivial contributions of the 
nucleon Dirac sea. 

In this paper we change the starting ideas of the RFMD having been 
formulated in [1,2]. We show that the RFMD is fully motivated by QCD 
and describes low-energy nuclear interactions in the nuclear phase of QCD 
through one-nucleon loop exchanges. Within the large Nc expansion [18,19] 
we justify the dominance of one-nucleon loop anomaly contributions to effec- 
tive Lagrangians describing the deuteron itself and processes of low-energy 
interactions of the deuteron coupled to other particles. 

^In Ref. [5] we have considered a modified version of the RFMD which is not well defined 
due to a violation of Lorentz invariance of the effective four-nucleon interaction describing 
N + N — > N + N transitions. This violation has turned out to be incompatible with 
a dominance of one-nucleon loop anomalies which are Lorentz covariant. Thereby, the 
astrophysical factor for the solar proton burning calculated in Ref. [5] and enhanced by a 
factor of 1.4 with respect to the recommended value [7] is not good established. 
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The paper is organized as follows. In Sect. 2 we discuss non-perturbative 
phases of QCD and formulate the RFMD from the first principles of QCD. 
In Sect. 3 we derive the effective Lagrangian for the deuteron field induced 
in the nuclear phase of QCD as the neutron-proton collective excitation 
(the Cooper np-pair) by a phenomenological local four-nucleon interaction. 
We demonstrate the dominant role of one-nucleon loop anomalies for the 
formation of the effective Lagrangian of the physical deuteron field. In Sect. 4 
we investigate the electromagnetic properties of the deuteron and derive the 
effective Lagrangian of the deuteron field, the Cooper np-pair, coupled to 
an external electromagnetic field describing the magnetic dipole and electric 
quadrupole moments of the physical deuteron. In the Conclusion we discuss 
the obtained results. 

2 Non— pert ur bat ive phases of QCD 

The derivation of the RFMD from the first principles of QCD goes through 
three non-perturbative phases of the quark-gluon system. We call them as: 
1) the low-energy quark-gluon phase (low-energy QCD), 2) the hadronic 
phase and 3) the nuclear phase. 

The low-energy quark-gluon phase of QCD can be obtained by integrat- 
ing over fluctuations of quark and gluon fields at energies above the scale of 
spontaneous breaking of chiral symmetry (SB^S) ~ 1 GeV. This results 
in an effective field theory, low-energy QCD, describing strong low-energy in- 
teractions of quarks and gluons. The low-energy quark-gluon phase of QCD 
characterizes itself by the appearance of low-energy gluon-field configura- 
tions leading to electric colour fluxes responsible for formation of a linearly 
rising interquark potential. The former realizes quark confinement. 

For the transition to the hadronic phase of QCD one should, first, in- 
tegrate out low-energy gluon degrees of freedom. Integrating over gluon 
degrees of freedom fiuctuating around low-energy gluon-field configurations 
responsible for formation of a linearly rising interquark potential one arrives 
at an effective field theory containing only quark {q) and anti-quark {q) de- 
grees of freedom. This effective field theory describes strong interactions at 
energies below the SB^S scale A^^ ~ 1 GeV. The resultant quark system pos- 
sesses both a chirally invariant and a chirally broken phase. In the chirally 
invariant phase the effective Lagrangian of the quark system is invariant un- 
der chiral U{3)xU{3) group. The chirally invariant phase of the quark system 
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is unstable and the transition to the chirahy broken phase is advantageous. 
The chirally broken phase characterizes itself by three non-perturbative phe- 
nomena: SB^S, hadronization (creation of bound quark states with quantum 
numbers of mesons qq, qqqq, baryons qqq and so on) and confinement. The 
transition to the chirally broken phase caused by SB^S accompanies itself 
with hadronization. Due to quark confinement all observed bound quark 
states should be colourless. As in such an effective field theory gluon degrees 
of freedom are integrated out, the entire variety of strong low-energy inter- 
actions of hadrons at energies below the SBxS scale ~ 1 GeV is described 
by quark-loop exchanges. 

Since nowadays in continuum space-time formulation of QCD the inte- 
gration over low-energy gluon-field configurations can be hardly performed 
explicitly, phenomenological approximations of this integration represented 
by different effective quark models with chiral U{3) x U (3) symmetry moti- 
vated by QCD are welcomed. 

The most interesting effective quark model allowing to describe analyti- 
cally both SBxS and bosonization (creation of bound qq states with quantum 
numbers of observed low-lying mesons) is the extended Nambu-Jona-Lasinio 
(ENJL) model [20-24] with linear [20,22] and non-linear [21,23] realization 
of chiral U{3) x U (3) symmetry. As has been shown in Ref. [24] the ENJL 
model is fully motivated by low-energy QCD with a linearly rising interquark 
potential and Nq quark colour degrees of freedom at Nc oo. In the ENJL 
model mesons are described as qq collective excitations (the Cooper gg-pairs) 
induced by phenomenological local four-quark interactions. Through one- 
constituent quark-loop exchanges the Copper gg-pairs acquire the properties 
of the observed low-lying mesons such as 7r(140), K(A90), 77(550), p(770), 
a;(780), i^*(890) and so on. For the description of low-lying octet and decu- 
plet of baryons the ENJL model has been extended by the inclusion of local 
six-quark interactions responsible for creation of baryons as qqq collective 
excitations [25]. 

Integrating then out low-energy quark-field fluctuations, that can be 

performed in terms of constituent quark-loop exchanges, one arrives at the 
hadronic phase of QCD containing only local meson and baryon fields. The 
couplings of low-lying mesons and baryons are described by Effective Chiral 
Lagrangians with chiral t/(3) x t/(3) symmetry [20-27]. 

The nuclear phase of QCD characterizes itself by the appearance of bound 
nucleon states - nuclei. In order to arrive at the nuclear phase of QCD we sug- 
gest to start with the hadronic phase of QCD and integrate out heavy hadron 
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degrees of freedom, i.e. all heavy baryon degrees of freedom with masses heav- 
ier than masses of low-lying octet and dccuplct of baryons and heavy meson 
degrees of freedom with masses heavier than the SB^S scale ~ 1 GeV. 
At low energies the result of the integration over these heavy hadron de- 
grees of freedom can be represented in the form of phenomenological local 
many-nucleon interactions. Following the scenario of the hadronic phase of 
QCD, where hadrons are produced by phenomenological local many-quark 
interactions as many-quark collective excitations, one can assume that some 
of these many-nucleon interactions are responsible for creation of many- 
nucleon collective excitations. These excitations acquire the properties of 
observed bound nucleon states - nuclei through nucleon-loop and low-lying 
meson exchanges. This results in an effective field theory describing nuclei 
and their low-energy interactions in analogy with Effective Chiral Lagrangian 
approaches [26,27]. Note that Chiral perturbation theory can be naturally 
incorporated into this effective field theory of low-energy interactions of nu- 
clei. 

We would like to emphasize that in this scenario of the quantum field 
theoretic formation of nuclei and their low-energy interactions nuclei are 
considered as elementary particles. For the first time, the representation of 
nuclei as elementary particles has been suggested by Sakita and Goebal [28] 
and Kim and Primakoff [29] for the description of electromagnetic and weak 
nuclear processes. We develop the quantum field theoretic approach to the 
interpretation of nuclei as elementary particles by starting with QCD. 

Following this scenario the deuteron, being the lightest bound nucleon 
state, appears in the nuclear phase of QCD as the neutron-proton collective 
excitation (the Cooper np-pair) induced by a phenomenological local four- 
nucleon interaction. Through one-nucleon loop exchanges [1-5,20-24] the 
Cooper np-pair with quantum numbers of the physical deuteron acquires the 
properties of the physical deuteron (i) the binding energy ed = 2.225 MeV, 
(ii) the electric quadrupole moment Qd = 0.286 fm^ [30] and so on. 

3 The deuteron as a Cooper np— pair 

In order to describe the deuteron induced as the Cooper np-pair we 
introduce a phenomenological local four-nucleon interaction caused by the 
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integration over heavy hadron degrees of freedom 



Ant(x) = -^JJ(X)/(X), (3.1) 

where gy is the phenomenological couphng constant of the RFMD [1-4], 
Mn = 940 MeV is the nucleon mass. We neglect here the electromagnetic 
mass difference for the neutron and the proton. As has been found in [1,2] 
the coupling constant gy is related to the electric quadrupole moment of the 
deuteron Q^: g^ = 2TT^Qr>M^ [2]. 

The baryon current j'^(x) with the quantum numbers of the deuteron is 
defined by [1-4] 

= -i []?(x)7''ra(a;) - ra=(x)7Xx)]. (3.2) 

Here p{x) and n{x) are the interpolating fields of the proton and the neutron, 
N'^{x) = C N'^{x) and N'^{x) = N'^{x)C, where C is a charge conjugation 
matrix and T is a transposition. In terms of the electric quadrupole moment 
of the deuteron the phenomenological local four-nucleon interaction Eq. ( |3.1|) 
reads 

C;M = -l^'QDjlix)fix). (3.3) 

Now let us discuss the behaviour of the phenomenological coupling constant 
gy/AM^ from the point of view of the large Nc expansion in QCD with the 
SU{Nc) gauge group at Nc — > oo [18,19]. Suppose, for simphcity, that the 
phenomenological four-nucleon interaction Eq.( p.l| ) is caused by exchanges 
of the scalar /o(980) and ao(980) mesons being the lightest states among 
heavy hadrons we have integrated out. 

Through a linear realization of chiral U{3) x f/(3) symmetry and the 
Goldberger-Treiman relation one can find that the coupling constant of cr- 
mesons (^o-nn, the gg-scalar mesons, coupled to the octet of low-lying baryons 
should be of order ^fg-NN ~ 0{y/Nc) at Nc oo- The scalar mesons /o(980) 
and ao(980) are most likely four-quark states with qqqq quark structure 
[31,32]. In the limit Nc — >■ oo such qqqq states are suppressed by a factor 
1/Nc [19]. Thus, an effective coupling constant of low-energy NN interaction 
caused by the qqqq scalar meson exchanges should be of order 0{1/Nc) at 
Nc —>■ oo. By taking into account that in QCD with Nc oo the nucleon 
mass Mn is proportional to Nc [19], Mn = NcMg, where Mg ~ 300 MeV is 
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the constituent quark mass, we can introduce the nucleon mass Mn in the 
effective couphng constant as a dimensional parameter absorbing the factor 
iV^, i.e. (7y/4M^. This is also required by the correct dependence of the 
deuteron mass on Nc- As a result the phenomeno logical coupling constant 
gy turns out to be of order 0{y^Nc) at Nc — > oo. 

We should emphasize that one does not need to know too much about 
quark structure of heavy hadron degrees of freedom we have integrated out. 
Without loss of generality one can argue that among the multitude of con- 
tributions caused by the integration over heavy hadron degrees of freedom 
one can always find the required local four-nucleon interaction the effective 
coupling constant of which behaves like 0{1/Nc) at Nc oo. As we show 
below this behaviour of the coupling constant of the phenomenological four- 
nucleon interaction is consistent with the large Nc dependence of low-energy 
parameters of the physical deuteron. 

The effective Lagrangian of the np-system unstable under creation of the 
Cooper np-pair with quantum numbers of the deuteron is then defined 

C'Pix) = n{x) {iYd^, - Mn) n{x) + p{x) {i^df, - Mn) p{x) 
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jt(a;)/(a;) + ..., (3.4) 



4M2 

where ellipses stand for low-energy interactions the neutron and the proton 
with other fields. 

In order to introduce the interpolating local deuteron field we should 
linearalize the Lagrangian Eq.( |3.4| ). Following the procedure described in 
Refs. [20-23] for the inclusion of local interpolating meson fields in the ENJL 
model we get 

^'^{x) n{x) {iYdf, - Mn) n{x) + p{x) {i^d^, - M^) p{x) 
+Ml Dl{x)D^{x) + gyjl{x)D^\x) + gynx)Dl{x) + . . . , (3.5) 

where Mq = 2Mn and D^{x) is a local interpolating field with quantum 
numbers of the deuteron. 

In order to derive the effective Lagrangian of the physical deuteron field 
we should integrate over nucleon fields in the one-nucleon loop approxima- 
tion [1-5,20-24]. The one-nucleon loop approximation of low-energy nuclear 
forces allows (i) to transfer nucleon flavours from an initial to a final nuclear 
state by a minimal way and (ii) to take into account contributions of nucleon- 
loop anomalies [1-5,17, 33-35], which are fully defined by one-nucleon loop 
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diagrams [33-35]. It is well-known that quark-loop anomalies play an im- 
portant role for the correct description of strong low-energy interactions of 
low-lying hadrons [20-27]. We argue the dominant role of nucleon-loop 
anomalies for the correct description of low-energy nuclear forces in nuclear 
physics. We demonstrate below the dominant role of nucleon-loop anomalies 
by example of the evaluation of the effective Lagrangian of the free deuteron 
field. 

The effective Lagrangian of the free deuteron field evaluated in the one- 
nucleon loop approximation is defined by [1-4]: 

J dt^xC^s{x) = j d^xM^Dl{x)D''{x) 

-jd^xj ^^^^ e-'k, . {X - X,) Dl{x)D^{x,) ^ W^k,; Q), (3.6) 

where the structure function U'^" {ki, Q) is given by 

r d^k f 1 1 1 

^''"{h-Q) = / — tr . . ^7^^ . .7^^ . (3.7) 

' J T^h [M^-k-Q-ki M^-k-Q I 

The 4-momentum Q — aki is a,n arbitrary shift of momenta of virtual nu- 
cleon fields with an arbitrary parameter a [1,2]. For the evaluation of the 
Q-dependence of the structure function n^^(A:i;(5) we apply the procedure 
suggested by Gertsein and Jackiw [32] (see also [2]): 

} d 

W'^iki; Q) - W'{ki- J dx — U^'ik,; xQ) = 



^ dx ^Q^T^M ^ ^7^^ -7 \ = 

i J T:h dk^ \M^-k-xQ-ki - k - xQ ] 

^2jdx ^lim /^^tT{{MN + k + xQ + kh^{M^ + k + xQ)Y]) = 

^ ^ ' 
= 2 {2Q^'Q^ - g2 g^"") + 2(A;f + k^Q^^ - ki ■ Q g'"'') = 
= -2 a(a + 1) (A;^ gi''' - 2 k^k^) . (3.8) 

Thus, we obtain 

W^{kr, Q) - W{kr, 0) = -2 a(a + 1) {kl gi"" - 2 k'^k\). (3.9) 
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We would like to emphasize that the r.h.s. of Eq.( ^.9| ) is an explicit expression 
completely defined by high-energy (short-distance) NN fiuctuations, since 
the virtual momentum k is taken at the limit k oo, and related to the 
anomaly of the one-nucleon loop diagram with two vector vertices (the VV- 
diagram) [33,35]. 

The structure function Il'^'^{ki] 0) has been evaluated in Refs. [1,2,5] and 
reads 

W{k,; 0) = ^{k^g^'' - k'^k-OMM^) + 2g^^[J,{M^) + MlUM^)], (3.10) 

where Ji(Mn) and J2{M^) are the quadratically and logarithmically diver- 
gent integrals [1,2,5]: 

r ct^k 1 d\k\k'^ 

r d^k 1 P d\k\k'^ , ^ 

The cut-off Ad restricts from above 3-momenta of low-energy fluctuations of 
virtual neutron and proton fields forming the physical deuteron [1-5]. Since 
in the RFMD the cut-off is much less than the mass of the nucleon, 
Mn ^ Ad [1,2], we use below the relation between the divergent integrals: 

Ji(Mn) = 2 J2(Mn) = t^- Oil/Nc). (3.12) 



Note that in Eg. ( p.lOD we have taken into account only the leading terms in 



the external momentum expansion, i.e. the fci-expansion. 

The justification of the dominance of the leading order contributions in 
expansion in powers of external momenta can be provided in the large Nq 
expansion. Indeed, in QCD with the SU{Nc) gauge group at A^^^ — > oo the 
baryon mass is proportional to the number of quark colours [19]: Mn ~ Nc- 
Since for the derivation of effective Lagrangians describing the deuteron and 
amplitudes of processes of low-energy interactions of the deuteron coupled to 
other particles all external momenta of interacting particles should be kept 
off-mass shell, the masses of virtual nucleon fields taken at Nc —>■ oo are 
larger compared with external momenta. By expanding one-nucleon loop di- 
agrams in powers of 1/Mn we get an expansion in powers of I/Nq- Keeping 
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the leading order in the large Nq expansion we are leaving with the leading 
order contributions in an external momentum expansion. We should empha- 
size that anomalous contributions of one-nucleon loop diagrams are defined 
by the least powers of an external momentum expansion. This implies that 
in the RFMD effective Lagrangians of low-energy interactions are completely 
determined by contributions of one-nucleon loop anomalies. The divergent 
contributions having the same order in momentum expansion are negligible 
small compared with the anomalous ones due to the inequality Mn ^ Ad 
[1-5] and the limit Nc — >■ oo. This justifies the application of the approx- 
imation by the leading powers in an external momentum expansion to the 
evaluation of the effective Lagrangians of the deuteron coupled to itself and 
an external electromagnetic field [2] , and the effective Lagrangians describing 
amplitudes of low-energy nuclear processes like the solar proton burning p 
+ p ^ D + e+ + i/e and so [2,4]. 

Collecting all pieces we get the structure function Ii^^{ki] Q) in the form 

U^^ih-Q) = ^(fciV'^-A;^^A;nJ2(MN)+2(7^^[Ji(MN)+M2j2(MN)] 

-2a{a + l){klg^'' -2k'ik^^). (3.13) 
The effective Lagrangian of the free deuteron field is then defined 



+ [m'o - ^ IMM^) + M^J,{M^)]^ Dl{x)D^{x), (3.14) 

where D^^{x) = d^D'^{x) — d^D^{x). We have dropped some contributions 
proportional to the total divergence of the deuteron field, since d^D^^^x) = 0. 
For the derivation of Eq. ( |3.14| ) we have used the relation 



/ ■ " "^'^ Dl{x)D.{x,){klg^'^ - k^.k'i) = 

= \j d'xDl,{x)D^^{x). (3.15) 

In order to get a correct kinetic term of the free deuteron field in the effective 
Lagrangian Eq.( |3.14| ) we should set [2] 



|^a(a + l) = l. (3.16) 
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Since a is an arbitrary real parameter, the relation Eg. ( p. 16 ) is valid in the 



case of the existence of real roots. For the existence of real roots of Eg. ( p. 161) 
the coupling constant gy should obey the constraint gy > Svr^ that is satisfied 
by the numerical value gy = 11.319 calculated at Nc = 3 [2]. Since gy ~ 
0{^/Nc) at Nc oo, Eg.( |3.16|) has real solutions for any Nc > 3. 



Due to Eg.( p.l6| ) the effective Lagrangian of the free deuteron field takes 
the form 




C^^ix) = -- f 1 + J2(Mn) 1 Dl,{x)D 



flU, 



-.2 



+ \M^ - ^ [MMn) + M^J2(Mn)] j Dlix)D'^ix). (3.17) 

By performing the renormalization of the wave function of the deuteron field 
[1,2] 

(^1 + ^J2(Mn)^ D^^ix) ^ D^^ix) (3.18) 

and taking into account that ^ we arrive at the effective Lagrangian 
of the free physical deuteron field 

£eff(x) = -^Dl{x)D^'^ix) + M^Dlix)D^ix), (3.19) 

where Md = Mq — £d is the mass of the physical deuteron field. The binding 
energy of the deuteron eb reads 



17 gl Ji(A/n) 17 



We have used here the relation between divergent integrals Eg.( |3.11 ) and 



expressed the phenomenological coupling constant gy in terms of the electric 
guadrupole moment of the deuteron g^ = 27i'^Qy)M^. 

At Nc oc the binding energy of the deuteron behaves like 0{1/Nc) 
as well as the electric guadrupole moment Qd and the coupling constant of 
the phenomenological local four-nucleon interaction Eg.( p.l|) . This testifies 
a self-consistency of our approach. Really, all parameters of the physical 
deuteron field are of the same order according to the large Nc expansion. 
This means that the vanishing of the coupling constant of the phenomeno- 
logical four-nucleon interaction Eg.( |3.1| ) in the limit Nc oo entails the 
vanishing of all low-energy parameters of the physical deuteron. 
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4 Electromagnetic properties of the deuteron 



The description of the deuteron as a Cooper np-pair changes a httle bit the 
description of the electromagnetic parameters of the deuteron given in Ref. [2]. 
We do not have more a "bare" deuteron field having magnetic dipole and elec- 
tric quadrupole moments. Therefore, both the magnetic dipole and electric 
quadrupole moments have to be described by the one-nucleon loop contri- 
butions. For the self-consistent description of the electromagnetic properties 
of the deuteron we cannot deal with only the baryon current j^{x) given by 
Eq.( |3.2|) and have to introduce the tensor current [2] 

J^''{x) = p^{x)af"'n{x) - n^{x)a'"'p{x), (4.1) 

where a''" = (YY - YY)/'^- 

The local four-nucleon interaction producing the deuteron as a Cooper 
np-pair reads now 



A„t(x) = -^jt(x)J'^(a;). (4.2) 



The baryon current J^{x) is defined by 

J^{x) = —i gy [p'^{x)Y^{x) ~ ^'^{x)l^p{x)] 

-7r^9,[P'i^)^''^^i^) - n-{x)cr''^p{x)], (4.3) 

where qt is a dimensionless phenomenological coupling constant [2]. The 
contribution of the tensor nucleon current looks like the next-to-leading term 
in the long-wavelength expansion^ of an effective low-energy four-nucleon 
interaction. 

The effective Lagrangian of the np-system unstable under creation of the 
Cooper np-pair with quantum numbers of the deuteron is then defined 

£°P(a;) = n{x){iYdf,- M^)n{x) +p{x){iYdf,- M^)p{x) 

jUx)nx). (4.4) 



4M2 



The linearalized version of the effective Lagrangian Eq. ( [4.4| ) containing the 
interpolating local deuteron field reads 

/:°P(x) ^ n{x) {lYd^ - Mn) n{x) + p{x) [lYd^, - Mn) p{x) 



^Due to proportionality Mn ^ Nc this expansion is related to the large Nc expansion. 
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+Ml DUx)D^{x) + g^3l{x)D^{x) + g^3^{x)Dl{x) 



+§-/,A^)D"'{^) + §-^J"'i^)DUi^), (4.5) 



where Mq = 2 Mn, D^{x) is a local interpolating field with quantum numbers 
of the deuteron and D^'''{x) = 9^L)^(x) - d''D''{x). 

The interactions with the tensor current give the contributions only to the 
divergent part of the effective Lagrangian of the free deuteron field defined 
now by [2] defined 



Aff(-) = ( - ^ + 1) + ^^ + '^3;? + '"^ J.(M,)) Dlix)D^^ix) 

+ (^0 - ^ IMM^) + M^MM^)]j Dl{x)D^{x). (4.6) 

Due to the relation Eq. (|3.16| ) the effective Lagrangian of the free deuteron 
field Eq. ([4.6|) takes the form 

+ [mS - ^ IMMn) + Mlh[M«)]j B* WB''(rc). (4.7) 

After the renormalization of the wave function of the deuteron field we arrived 
at the effective Lagrangian defined by Eq.( |3.19| ) with the binding energy of 
the deuteron depending on gy and gT [2] 

17 g^ Ji(Mn) / 48 2A g^ _ 
^° 48 772 Mn 17 gy 17 g'vj 

^ i^g.A^(l + ^-+^4), (4.8) 

18 V 17 gy 17 g^ J 

where we have used the relation between divergent integrals Eq.( |3.11| ) and 
expressed the phenomenological coupling constant gy in terms of the electric 
quadrupole moment of the deuteron g^ = 2'7i'^Q]:,M^. In order to make 
the prediction for the binding energy much more definite we have to know 
the relation between the phenomenological coupling constants gy and gT- 
For this aim we suggest to consider the electromagnetic properties of the 
deuteron. 
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Including the electromagnetic field by a minimal way d^ + ie A^{x) 
where e and A^{x) are the electric charge of the proton and the electromag- 
netic potential we bring up the linearalized version of the Lagrangian Eq, 
to the form 



n{x) {i^^'d^ - Mn) n{x) + p{x) [i^^'d^ - Mn) + Dl{x)D^'{x) 

Mo ^ ' Mo '^"^^'^ 

9t 



+9v3l{x)D^{x) + gY3'{x)Dl{x) + ^/,.D^'{^) + ^^"'^i. 



ep{x)^^p{x)A,{x)-te^JlAx){A^{x)D'^{x)-A'^{x)D^{x)) 
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+ie^.r^{x){A^{x)D,{x) - A-{x)D^{x)). (4.9) 

The effective Lagrangian describing both the magnetic dipole and electric 
quadrupole moments of the deuteron has been evaluated in Ref. [2] and reads 

5C^^^,{xU = ,e ^"^^~/^^ Dlix)A^ix)D^ix) 



— te- 



6 TT 

. 4:ag^-bg^ 



6 TT 
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^D^''{x)A,{x)Dl{x) 



+ ze^{2b + 3)Dl{x)D,{x)F^''{x) 

+ ^e{l + a)^^Dl{x)D^\x)F,'^ix), (4.10) 

where F^^{x) = dfj,A,y{x) — di^A^^x) is the electromagnetic field strength, a 
and b are arbitrary parameters related to the ambiguities of the one-nucleon 
loop diagrams with respect to a shift of virtual nucleon momentum. We 
consider them as free parameters of the approach [2]. 

In order to fix parameters it is convenient to write down the total effective 
Lagrangian of the physical deuteron coupled to an external electromagnetic 
field 



>CgLM(^)eff = -lDl,ix)D^'^ix) + M^Dlix)D^ix) 
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(4.11) 



Two terms having the structure D'^'^ {x) A^{x) Dl{x) and Djjy(x)y4'*(a;)D'^(a;) 
should describe the interaction of the deuteron with an external electromag- 
netic field included by a minimal way, whilst the last two terms are respon- 
sible for the non-trivial contributions to the magnetic dipole and electric 
quadrupole moments of the deuteron. In terms of the parameters of the ef- 
fective interactions Eg. ( [4.11D the magnetic dipole moment /id, measured in 
nuclear magnetons, and the electric quadrupole moment Qd, measure in fm^, 
of the deuteron are given by 



Qd 



+ a + b) 



at the constraint 
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(4.13) 



reducing the first two terms in effective Lagrangian Eq. (f4.10 ) to the standard 
minimal form which can be obtained from the effective Lagrangian of the free 
deuteron field by the shift d^D^{x) {d^ + i e A^[x))Dy{x). 

By retaining the former relation between the electric quadrupole moment 
and the coupling constant gy, Qb = 2g'y/7i'^M^ [2], one can show that the 
experimental values of the magnetic dipole moment of the deuteron /id = 
0.857 and the electric quadrupole moment Qn = 0.286 fm^ can be fitted by 
the following values of the parameters a and b: a = — 0.442 and b = — 4.418 
at gy = 11.319 [2]. This gives the relation between the coupling constants 
gy and gT 



gT = - 1.662 5(v. 



(4.14) 



Substituting this relation into Eq.( [1.8|) we can describe the experimental 
value of the binding energy of the deuteron = 2.225 MeV at the cut-off 
Ad = 115.729 MeV. The spatial region of virtual nucleon field fiuctuations 
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forming the physical deuteron related to this value of the cut-off 1/Ad ~ 
Pd = 1.705 fm agrees good with the effective range of the deuteron pd = 
(1.759 ± 0.005) fm [6]. 

The effective Lagrangian of the deuteron field coupled to an external 
electromagnetic field is given by 

^^LM{x)eS = — [(d^ - ieA^(x))Di{x) - (d, - ieA,{x))Dl(x)] 

X [{df'+ ieA^'{x))D''{x) - {d''+ ieA''{x))D''{x)] + M^Dl{x)D^'{x) 

- i e /XD Dlix) D,{x) F^''{x) + teQ^ Dl{x) D''\x) F/(a;), (4.15) 

where Ji^ = 14.682 and Qd = 0.514 Qd. The term of order O(e^) can 
be also derived in the RFMD by using shift ambiguities of one-nucleon loop 
diagrams. This term is required by the electromagnetic gauge invariance of 
the effective Lagrangian of the deuteron field, but it does not affect on the 
electromagnetic parameters of the deuteron which are of order 0{e). 

Conclusion 

Unlike the first formulation of the RFMD given in Refs. [1,2], where 
we have stated that the RFMD is far from being induced by the dynamics 
of QCD and seems like an old-fashion approach [2], in this paper we have 
shown that the RFMD is in complete agreement with QCD and can be for- 
mulated from the first principles of QCD. The RFMD describes low-energy 
nuclear forces in the nuclear phase of QCD in terms of one-nucleon loop ex- 
changes. One-nucleon loop exchanges provide a minimal way of the transfer 
of nucleon flavours from an initial to a flnal nuclear state and allow to take 
into account contributions of nucleon-loop anomalies. These anomalies are 
related to high-energy fluctuations of virtual nucleon flelds, i.e. the NN 
fluctuations, and fully determined by one-nucleon loop diagrams [33-35]. 
The dominance of contributions of one-nucleon loop anomalies to effective 
Lagrangians describing low-energy interactions of the deuteron coupled to 
itself and other particles is justifled in the large Nc expansion in QCD at 
Nc oo. It is well-known that anomalies of quark-loop diagrams play an 
important role for the correct description of strong low-energy interactions 
of low-lying hadrons. We argue an important role of nucleon-loop anomalies 
for the correct description of low-energy nuclear forces in the nuclear physics. 
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It is interesting that nucleon-loop anomalies can be interpreted as non- 
trivial contributions of the non-perturbative quantum vacuum - the nucleon 
Dirac sea [17]. In nuclear physics the influence of the nucleon Dirac sea for 
low-energy properties of finite nuclei has been analysed within quantum field 
theoretic approaches in the one-nucleon loop approximation [16]. Unfortu- 
nately, in these approaches contributions of one-nucleon loop anomalies have 
been taken into account. The RFMD allows to fill this blank. 

For the derivation of the RFMD from the first principles of QCD we distin- 
guish three non-perturbative phases of QCD: 1) the low-energy quark-gluon 
phase (low-energy QCD), 2) the hadronic phase and 3) the nuclear phase. 
Skipping over the intermediate low-energy quark-gluon phase by means of 
the integration over high- and low-energy quark and gluon fiuctuations one 
arrives at the hadronic phase of QCD containing only local hadron fields 
with quantum numbers of mesons and baryons coupled at energies below the 
SBxS scale ~ 1 GeV. The couphngs of low-lying mesons with masses 
less than the SB^S scale to low-lying octet and decuplet of baryons can be 
described by Effective Chiral Lagrangians with chiral U{3) x U{3) symmetry. 

Integrating in the hadronic phase of QCD over heavy hadron degrees of 
freedom with masses exceeding the SBxS scale one arrives at the nuclear 
phase of QCD which characterizes itself by the appearance of bound nucleon 
states - nuclei. At low energies the result of integration over heavy hadron 
degrees of freedom can be represented in the form of phenomenological local 
many-nucleon interactions. Some of these interactions are responsible for 
creation of many-nucleon collective excitations which acquire the properties 
of observed nuclei through nucleon-loop and low-lying meson exchanges. 
This effective field theory describes nuclei and processes of their low-energy 
interactions by considering nuclei as elementary particles [28,29]. 

Following this scenario of the description of nuclei and their low-energy 
interactions from the first principles of QCD the deuteron should be produced 
in the nuclear phase of QCD by a phenomenological local four-nucleon in- 
teraction as the Cooper np-pair with quantum numbers of the deuteron. 
The properties of the physical deuteron, i.e. the binding energy, the elec- 
tric quadrupole moment and so, the Cooper np-pair acquires through one- 
nucleon loop exchanges. As has been shown the main part of the kinetic 
term of the effective Lagrangian of the free physical deuteron field is induced 
by the contribution of high-energy (short-distance) fluctuations of virtual 
nucleon flelds related to the anomaly of the one-nucleon loop diagram with 
two vector vertices - the VV-diagram. 
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In turn, the magnetic dipole /iD and electric quadrupole moments of the 
physical deuteron Qd are fully determined by high-energy (short-distance) 
fluctuations of virtual nucleon fields related to the anomalous contributions of 
the one-nucleon loop diagrams with three vector vertices (the 7yy-diagram) 
[2]. Thus, high-energy (short-distance) fluctuations of nucleon fields related 
to anomalies of one-nucleon loop diagrams play a dominant role for the 
formation of the physical deuteron from the Cooper np-pair. 

As regards low-energy (long-distance) fluctuations of virtual nucleon 
fields they give a significant contribution only to the binding energy of the 
deuteron ^d- The strength of low-energy (long-distance) fiuctuations of vir- 
tual nucleon fields is restricted by the cut-off Ad = 115.729 MeV. The spatial 
region of virtual nucleon field fiuctuations forming the physical deuteron re- 
lated to this value of the cut-off 1/Ad ~ Pd = 1.705 fm agrees good with the 
effective range of the deuteron pD = (1.759 ± 0.005) fm [6]. 

It is well-known that in the potential model approach to the description 
of the deuteron the electric quadrupole moment of the deuteron is caused 
by nuclear tensor forces which play an important role for the existence of the 
deuteron as a bound np-state. 

The proportionality of the coupling constant of the phenomenological 
local four-nucleon interaction Eq. (p.3|) , responsible for creation of the Cooper 
np-pair with quantum numbers of the deuteron, and the binding energy of 
the deuteron £d Eq.( p.20|) to the electric quadrupole moment Qd testifies an 
important role of nuclear tensor forces for the formation of the deuteron in 
the RFMD. 

To the evaluation of one-nucleon loop diagrams defining effective La- 
grangians describing processes of low-energy interactions of the deuteron 
coupled to itself and other particles we apply expansions in powers of the 
momenta of interacting particles and keep only leading terms of the expan- 
sions [1-5]. This approximation can be justified in the large Nc expansion. 
Indeed, in QCD with the SU{Nc) gauge group at Nc — > oo the nucleon mass 
is proportional to the number of quark colours [19]: Mn ~ Nc- Since for the 
derivation of effective Lagrangians describing the deuteron and amplitudes 
of low-energy nuclear processes all external momenta of interacting particles 
should be kept off-mass shell, the masses of virtual nucleon fields are larger 
compared with the external momenta. An expansion of one-nucleon loop 
diagrams in powers of 1/Mn giving an external momentum expansion corre- 
sponds to the expansion in powers of I/Nq- In this case the leading order in 
the large Nq expansion gives the leading order contributions in the expansion 
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in powers of external momenta of interacting particles. We should emphasize 
that anomalous contributions of one-nucleon loop diagrams are determined 
by the least powers in external momentum expansions. Thereby, the dom- 
inance of contributions of nucleon-loop anomalies to effective Lagrangians 
describing low-energy nuclear processes in the RFMD is fully supported by 
the large Nc expansion. The accuracy of this approximation is rather high. 
Indeed, the real parameter of the expansion of one-nucleon loop diagrams is 
1/M^ ~ l/iV^ but not 1/Mn ~ 1/Nc- Thereby, next-to-leading corrections 
should be of order 0{1/N^). 

The inclusion of the interaction of the deuteron field with the tensor nu- 
cleon current Eq. (|4.1| ) has given a possibility of the self-consistent description 
of the electromagnetic properties of the deuteron, the magnetic dipole mo- 
ment fij) and the electric quadrupole moment Qd, in terms of effective inter- 
actions of the Corben-Schwinger and Aronson kinds induced by one-nucleon 
loop diagrams [2]. By fitting the experimental values of the magnetic dipole 
moment fi^ = 0.857yUN, where /iN = e/2M^ is a nuclear magneton, and the 
electric quadrupole moment = 0.286 fm^ supplemented by the require- 
ment of the electromagnetic gauge invariance of the effective Lagrangian of 
the deuteron field coupled to an external electromagnetic field [2] we have got 
the relation between the coupling constants and ^^t: = — 1.662 (yfy. Due 
to this relation the experimental value of the binding energy of the deuteron 
can be described by the cut-off Ad = 115.729 MeV. This corresponds to 
the spatial region of virtual nucleon field fluctuations forming the physical 
deuteron 1/Ad ~ Pd = 1.705 fm agreeing good with the effective range of 
the deuteron pD = (1.759 ± 0.005) fm [6]. As has been stated in Ref. [2] the 
contributions of the tensor nucleon current Eq.( ^.l| ) to the amplitudes of low- 
energy nuclear processes of astrophysical interest such as the neutron-proton 
radiative capture n -|- p D + 7 for thermal neutrons and the solar proton 
burning p-f-p— >-D-|-e^ + z/e can be neglected. The same statement is 
valid for the disintegration of the deuteron by anti-neutrinos and neutrinos 
induced by charged Pe + D^e^ + n + n, i/g + D— i>e^ + p + p and neutral 
i^ei^e) + D — > Uci^c) + n + p weak currents, and the pep process, p + e~ + 

p ^ D + I/e. 

The quantum field theoretic scenario to treating nuclei as many-nucleon 
collective excitations induced by phenomenological local many-nucleon in- 
teractions allows a plain extension of the RFMD by the inclusion of light 
nuclei ^He, and ^He as three- and four-nucleon collective excitations. 
The binding energies and other low-energy parameters of these excitations 
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should be defined through nucleon-loop and low-lying meson exchanges. The 
extension of the RFMD by the inclusion of '^He, "^H and '^He would give a 
possibility to continue investigations of the reactions of the p-p chain [36] 
started with the reaction p + p— >D + e"'" + t/e and to apply the extended 
version of the RFMD to the description of the reactions p + D — > ^He + 7, 
p + ^He — > ^He + e'^ + and so on. 

We would like to emphasize that Chiral pcrtTirbation theory can be nat- 
urally incorporated into this effective field theoretic approach to physics of 
low-energy interactions of nuclei in terms of Effective Chiral Lagrangians 
with chiral U{3) x U{3) symmetry describing low-lying baryons and mesons 
interacting at low energies [37]. Unfortunately, the discussion of the inclu- 
sion of Chiral perturbation theory into the RFMD and the description of 
low-energy nuclear processes in the RFMD goes beyond the scope of this 
paper. Therefore, we relegate readers to Refs. [37,38], where these problems 
and a comparison of the RFMD with other approaches have been discussed 
in details. 
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